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INTRODUCTION: Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, 
particularly in the lower extremities such as the proximal femur [1]. The most commonly used countermeasure 
against bone loss has been prescribed exercise [2]. However, current exercise countermeasures do not completely 
eliminate bone loss in long duration, 4 to 6 months, spaceflight [3,4], leaving the astronaut susceptible to early onset 
osteoporosis and a greater risk of fracture later in their lives. The introduction of the Advanced Resistive Exercise 
Device, coupled with improved nutrition, has further minimized the 4 to 6 month bone loss. But further work is 
needed to implement optimal exercise prescriptions [5]. In this light, NASA’s Digital Astronaut Project (DAP) is 
working with NASA physiologists to implement well-validated computational models that can help understand the 
mechanisms of bone demineralization in microgravity, and enhance exercise countermeasure development. 

METHODS: The objective of the DAP computational modeling effort is to enable simulations in time of changes in 
bone mineral density (BMD) and Bone Volume Fractions (BVF) under the conditions of skeletal unloading and 
changes in physiological processes encountered in microgravity. Since the geometry of the remodeling units or bone 
packets that are removed and replaced during remodeling differ, separate modules for trabecular bone and cortical 
bone are developed. Key elements of the computational model include: Bone resorption (formation) rate varies with 
activation density, volume of remodeling unit removed (replaced), and active osteoclast (osteoblast) population. The 
active osteoblast and osteoclast populations vary according to the cellular dynamics mediated by hormones, 
proteins, ligands and receptors. The well-known adaptive response theory of Frost drives the bone response to 
variations in skeletal loading. Within the expressions for rates of changes of the cellular populations, assumptions 
for the ligand and receptor expressions are modeled in accordance with the American Society of Bone and Mineral 
Research educational literature. 

INITIAL MODEL DEVELOPMENT The model’s initial development focuses on the femoral neck. Remodeling 
unit dimensions for the femoral neck (particular for cortical bone) were identified in the literature to make the model 
specifically applicable to the femoral neck, although many other model parameters were based on general bone 
knowledge. For model validation, we used BMD changes of control subjects in the current 70 bed rest study and 
available data from the 17-week bed rest studies conducted in the past (Figure 1). Volumetric bone densities for the 
70 bed rest were obtained at pre and post bed rest via Quantitative Computed Tomography (QCT). However, bone 
densities from past bed rest studies were obtained via Dual-energy X-ray Absorptiometry (DXA). Given that DXA 
is 2-D integrated cortical and trabecular, and the computational model tracks BVF changes, the DXA values were 
mapped to equivalent QCT integral volumetric density values in order to run simulations with the DXA data. Our 
poster will discuss in detail how the model tracks BVF and the preliminary model validation results. 



DISCUSSION AND FUTURE WORK: Our results show that a good foundation has been laid for establishing a 
physiologically accurate bone remodeling model. For example, mean BMD data for the 70-day bed rest is well 
within the 95% confidence interval of the model prediction. Future work will integrate the bone remodeling model 
with DAP’s biomechanical exercise models to predict the benefits of exercise induced load stimulus from different 
exercise prescriptions for maintaining bone at the femoral neck. The model will also be extended to include 
predictions for the lower lumbar spine, calcaneus, trochanter and the integrated proximal femur. 

REFERENCES: [1] Buckley (2006) Space Physiology. Oxford University Press, New York; [2] Layne and Forth 
(2008) Aviat. Space Environ. Med. 79, 787-794; [3] Lang et al. (2004) J. B Min Res. 19(6), 1006-1012; [4] Keyak et 
al. (2009) Bone 44 (2009) 449-453; [5] Smith et al. (2012) J. Bone & Min Res. 
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Background 


Bone Remodeling Model Implementation Plan 


General Description of the DAP Bone Remodeling Model 


Why Quantifying Change in Bone via Bone Remodeling is Objective of NASA Digital Astronaut Project (DAP) 

• One of the main objectives is to provide a tool to help HHC address Bone Gap Osteo 4 : We don’t know the contribution of 
each risk factor on bone loss and recovery of bone strength and which factors are the best targets for countermeasure 
application; and Osteo7 : We need to identify options for mitigation of early onset osteoporosis before, during, and after 
spaceflight. 

■ Skeletal loading along with endocrine regulation and local biochemical mediators are what drives the physiological 
mechanism of bone remodeling to maintain bone. 

■ Exercise induced loading, with appropriately input to a model can approximately predict the effect of specific exercise 
prescription and thus help to evaluate its benefits as a countermeasure option. Integrates with DAP Biomechanics 
Model and the DAP Muscle Model. 


Importance for the New Finite Element Based Strength Standard 


Other main objectives intend to inform the HHC Bone Discipline's efforts to address Bone Gap Fracture 3 . We need a 
validated method to estimate the Risk of Fracture by evaluating the ratio of applied loads to bone fracture loads for 
expected mechanically-loaded activities during and after a mission 

■ One effort is underway to evaluate Finite Element fFE) estimates of bone strength (aka bone fracture loads) as a 
potential standard for bone health. 

■ A bone remodeling formulation that quantifies dynamic changes in bone has the potential of tracking changes in 
volume fractions that can relate to QCT BMD and ash density estimates, upon which FE bone strength is based [1], In 
addition coupling a BR model with a QCT based FE model may also provide geometry changes. 



What does it do? 

It tracks changes in the bone when the balance between formation and 
resorption in the bone remodeling process becomes unbalanced. 


How does it do it? 

The cellular physiology, remodeling unit mechanisms, and mechano- 
transduction theory that drive the process are described mathematically. 


How does the computational algorithm work? 

Rates of change of bone volume fraction and cell populations are set to zero 
(Balanced healthy state with steady bone density). 


Balance is broken by skeletal unloading, and rate of change is no longer 0. 


The system including bone properties and cell populations are integrated in 
time to estimate the change. 


NOTE: Model parameters and methodology are currently focused on the femoral neck. 


Mathematical Description 


Modeling the Influence of Skeletal loading 


Model Representation of Bone Loss Due to Insufficient Mechanical 
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State Variables and Definitions 


Symbol Definitions in the Cell Equations 


via mass balance equations. The 

dynamics is a considerable 
modification of the work of 
Lemaire et al. [2] and Pivonka et 
al. [3] with the addition of 
effectors related to skeletal 


The most likely intermediaries that enable sensor cells to trigger effector cells is NO and PGE-2 151- 


Released by Osteocvtes and Osteoblasts under mechanical stimulation 

"Stimulates production of OPG 
Inhibits production of RANKL 


Mediates differentiation of osteoblasts 
induced by TGF - p 


Stimulates proliferation of osteoblasts 
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The model gages the level of expression of NO and PGE 2 according to the level of bone apposition o 
bone resorption suggested by the daily strain £ in Frost’s Mechanostat Theory as outlined below: 


Sensing strength or response level (SL) 
defined in relation to bone strain 


sl = /(*)= F? )Fr +i 


Complete Unloading £ = 0 SL = 0 
Remodeling Balance £ = s 0 SL = 1 
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Mathematical model of the Mechanostat. 


NOTE: Osteocytes are generally understood to be the sensor cells 

Rate per cell 


PGE Production Rate - 


Production rate of NO and 
PGE 2 per cell are defined to 

be proportional to SL NO Production Rate ■ 


+S a =p a xSLxY d xBVF 


*-S N = p N x SL x YJx BVF 

t 

Rate per cell 



Converting Experimental Data to Model Variables 


Preliminary Validation Results for Bone Deconditioning Simulations 


Future Work 


A Method for Mapping vBMD to BVF 


Apparent (dry) density p m 


Given: 

A Pre Bed Rest QCT BMD value. 

A Bed Rest Duration Length of N days. 
A Post Bed Rest QCT value. 


Ash fraction a = 


morgamc mass + organic mass 
D m - density Mineralizd bone D 0 - density of Osteoid 


Ash fraction has a theoretical limit 0.7 [4], 
T »nr_ apparentdensity 


true tis sue density p, 
M ■ D m ■ 0.7 M-D„ 
~ M-D m +0-D c OT M 


p,(glcc) = lAl 
+ OD a 


.29 a [4] 

= p, (Usedin Model) 


1. (a) Convert p QCT to p ash (e.g. Keyak regression) 

(b) Convert p mh to p app (e.g. Schileo regression) 

(c) Compute initial ash fraction a = p mh / p app 

2. Initial value M = p mh / (0.7 x D m ) 

Solve for initial value O using O. definition. 

3. Run computational simulation subject to loading 
history (i.e. bed rest) for N days to track change 
in M, O, a, p mh , p QCT (BMD), and BVF 

4. Compare BMD to QCT BMD 
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Preliminary validation results for predicting lc 

(a) Trabecular bone 

(b) Cortical bone After 70 days of bed re 


(c) Time course change of mean Dexa BMD for 18 control 
subjects during 17 weeks of bed rest. 


Near Term: 

• Develop/formulate a daily load formula for quantifying exercise 
induced loading and test against exercise treated subjects (e.g. 
CFT70 study) 


Long Term : 

• Develop method for transforming force data from biomechanics 
modeling of specific exercise devices into stress/strain input 

• Integrate the computational model with Finite Element Method 

• Validate model using QCT data from spaceflight research 

• Develop model for predicting bone adaptation for trochanter, total 
proximal femur and lower lumbar 

• Bone adaptation prediction for more than 180 days of spaceflight 
exposure with exercise countermeasure 
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